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Solutions to Activities 
Solutions to Exercises 


You should schedule four study sessions for your work on this 
chapter. The chapter includes a lot of computer work — more 
than any of the other chapters in this block. You will need to 
use a computer for three of the four sessions — the second, 
third and fourth. Each section of the chapter depends on 
ideas and results from the preceding sections or uses features 
of the software introduced in earlier sections, so we 
recommend that you study the sections in order. A possible 
study pattern is as follows. 


Study session 1: Sections 1 and 2. 


Study session 2: Section 3. You will need access to your 
computer for this session, together with the statistics software 
package, Computer Book D and an audio tape player. 


Study session 3: Section 4. You will need access to your 
computer for this session, together with the statistics software 
and Computer Book D. 


Study session 4: Sections 5 and 6. You will need access to 
your computer, together with the statistics software and 
Computer Book D for Section 5, 


If you follow this study pattern, then the second session will 
probably be the longest. 


Introduction 


Many commonplace phenomena exhibit variability — the lengths of 
pregnancies, the heights and weights of adults, the amounts people earn, 
annual rainfall, the times between major earthquakes or other natural 
disasters, and so on. The possibilities are too numerous to count. What 
they all have in common is that none can be predicted with certainty. 
For instance, when the next major earthquake will occur is unknown; 
sometimes two or three occur within a few months of each other, while at 
other times there will not be one for a year or more ~ the times between 
major earthquakes vary. And although human pregnancies are supposed 
to last 40 weeks, few babies arrive exactly ‘on time’ ~ the lengths of 
pregnancies vary. Babies’ birthweights vary too. And some men and 
women are taller than others. 


For each of these situations, we might collect some relevant data, calculate 
the mean, and use this to estimate, for instance, when the next major 
earthquake might occur, or the average height of all men, or the average 
weight of a baby at birth. But this will not tell us anything about how 
inaccurate our estimate might be, and it could be important to know this 
also. 


There are many situations where knowledge of how physical measurements 
vary is important. Clothing designers need to know about the variability of 
many measurements such as, for example, chest, bust, waist and hip girth, 
arm length and inside and outside leg length; and they also need to know 
how these measurements are related to each other. Manufacturers need to 
know the distributions of such measurements in order to make informed 
decisions about the numbers of each size to produce. 


Because humans vary in size, the distributions of physical measurements 
are also of great importance in the design of equipment for human use. 
The first systematic attempt to develop principles for equipment design 
and to collect data and apply the information obtained took place during 
the Second World War. This knowledge was used most effectively in the 
design of aircraft cockpits, oxygen masks and clothing. 


In the late 1940s, investigations were carried out in the United States into 
the use of all types of equipment, both civilian and military, with a view 
to improving safety. Among other things, data were collected on 
measurements such as sitting eye height, arm reach and leg length for bus 
and truck drivers. Poor design increases the risk of accidents: for example, 
arms tire if the steering wheel is too far from the body, and if the foot 
pedals are too close, then excessive bending of the ankles is required. 
Effective use of equipment depends on good design, and this requires 
knowledge of the physical dimensions of those likely to operate the 
equipment. 


Ideally, all men and women should be able to operate all equipment: car 
manufacturers, for instance, want their vehicles to be usable by everyone. 
And, in practice, provided that the overall dimensions are chosen correctly, 
adjustable devices can accommodate almost everyone. However, relevant 
data on the distribution of the physical dimensions of people are essential 
for the designers if the goal of universal operability is to be achieved. 
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The design of equipment such as aircraft cockpits or driving cabs or cars is 
very complicated, so, in order to introduce and develop ideas about how to 
model variation, we shall consider the much simpler problem of how to 
model the heights of men. 


This chapter is about modelling the variation observed in data; and it is 
concerned primarily with one particular model for variation called the 
normal distribution. In Section 1, we begin by discussing the general 
problem of how to model the variation observed in a sample of data, and 
then the normal distribution is introduced. In Chapter D1, the geometric 
distribution was introduced as a model for the variation observed in the 
number of trials of an experiment needed to obtain a success. You saw 
that, in general, choosing a geometric distribution to model a particular 
situation involved fixing the value of one parameter: p, the probability of 
success in each trial of the experiment. Choosing a normal distribution to 
model the variation observed in data involves fixing the values of two 
parameters. In Section 2, we discuss what these parameters represent and 
how data are used to select suitable values for them. The ideas introduced 
in the first two sections are used in Sections 3 and 4 to select models for 
the variation in a number of phenomena ~ in each case, the model is 
chosen to reflect the variation observed in a sample of data. 


There is considerable evidence that, within a country, human physical 
measurements change over a period of time. For example, records show 
that the average height of American soldiers was 0.7 inches greater in the 
Second World War than in the First World War. And in the 1960s, the 
trend for average heights to increase was noted for many other Western 
countries and also for Japan. In Section 4, you are invited to explore a 
large data set collected at the beginning of the 20th century on heights of 
fathers and sons, and to investigate whether, in general, the sons were 
taller than their fathers. Finally, the properties of normal distributions are 
investigated in Section 5 and summarised in Section 6. Sections 3, 4 and 5 
each require you to have access to your computer. 


There are two main learning file themes for this chapter. The first involves 
you in using the software package relevant to this block. This is the data 
analysis software OUStats for MS'T121, which from now on we shall call 
OUStats. An important part of this block is for you to use OUStats to help 
your learning and exploration of the statistical ideas introduced. But for 
you to be able to achieve this aim, you also need to become familiar with 
the software package. Many software packages have features in common, so 
skills developed and acquired with one package can often be transferred to 
use with another. Mathcad and OUStats are both Windows-based 
software packages, and some operations are carried out in a similar way for 
the two packages. Thinking about how you approach learning and using a 
new computer package can be very informative, so that you can use your 
experience from one situation to help you in a new one. 


The other theme continues to develop and consolidate your communication 
skills. In Chapter D1, you were encouraged to think about and practise 
discussing ideas — both from the viewpoint of how discussion work can be 
helpful for learning and because talking ideas through with people is often 
a necessary part of working in a team and completing a project. In this 
chapter, we take this one step forward and consider more formal 
presentations. Many areas of work and study involve communicating 
verbally in a more formal way by making a structured oral presentation. 


INTRODUCTION 


If you feel this has no relevance to you, take a few minutes to consider 
what this type of activity involves. To make an effective presentation, you 
need to have a clear, logical structure for the work you present, and an 
awareness of what the audience needs to know. It is also important to 
consider the effect of the talk on the audience’s understanding of the 
problem and its solution. In this way, the presenter is forced to look at the 
problem from different perspectives. Both preparing and giving a 
presentation can be very effective ways of assessing your own 
understanding of the work. 


Activity 0.1 Using experience for new learning 


As you begin to use OUStats, try to remember how you began using 
Mathcad. Is there anything you would do differently? What worked well, 
and were there aspects that worked less well as you became familiar with 
the package? What advice would you give to another student beginning 
work with Mathcad, in terms of how to get started and how to use the 
package in an effective way? Make notes on Learning File Sheet 1. 


Taking time to consider what worked well and which aspects of your work 
you wish to change can be a useful way to get involved in self-assessment: 
and so help you to improve your learning and performance. 


1 Models for variation 


In Chapter D1, you saw how to model the variability in one particular 
quantity — the number of rolls of a die needed to obtain a six. By assuming 
that the die is ‘fair’ — that is, that each score is equally likely to occur on 
each roll of the die — we developed a model for the variation in the number 
of rolls required to obtain a six: this model is the geometric distribution. 
The model specifies the probability of obtaining each possible number of 
rolls: 1,2,3,.... It gives a good indication of the likely variation in the 
number of rolls that may be needed to obtain a six on different occasions. 
This distribution is illustrated in Figure 1.1(a). 


In Section 2 of Chapter D1, you saw that simulations, based on the 
assumption that the six possible scores on each roll are equally likely to 
occur, gave rise to frequency diagrams which were similar in shape to 
the probability distribution in Figure 1.1(a). An example is shown in 
Figure 1.1(b), showing that the model seems to be a good one for what 
actually happens. 


Probability 4 
0.24 
O41 
10 20 30 
(a) Number of rolls needed for a six 
Frequency 4 


(b) Number of rolls needed for a six 


Figure 1.1 (a) A geometric distribution with p=} (b) The results of a 
simulation 


Now consider how we might model the variation in men’s heights. Since 
it would not be practicable to measure the heights of all men, the first 
step is to obtain some data. It proved surprisingly difficult to find data 
on the heights of individual men. Data on the distributions of human 
measurements can be of considerable commercial value, so many 
organisations that collect such data are reluctant to allow their 
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publication. However, some searching turned up an old data set which 
contains the heights of 1000 Cambridge men in 1902. These heights are a 
sample of the heights of all Cambridge men in 1902, and will be used to 
illustrate how the variation observed in men’s heights may be modelled. 


Figure 1.2 shows a frequency diagram for the heights of 1000 Cambridge 


men in 1902. The heights were recorded to the nearest inch; so, for In 1902, inches were used in 
instance, all men strictly between 69.5 and 70.5 inches tall were recorded the UK to measure heights, 
as 70 inches tall and are represented by the bar centred on 70 inches. rather than metric units. 


(Rounding of heights was used: a height measured as exactly 69.5 inches 
was recorded as 70 inches; a height measured as exactly 70.5 inches was 


recorded as 71 inches.) These data are from an article in which various Source: W. R. MacDonell, 
measurements taken on convicted criminals were compared with ‘On criminal anthropometry 
corresponding measurements for men in the general population. Around and the identification of 
this period, there was considerable interest in discovering whether there criminals’, Biometrika 1 


were differences between various physical measurements for criminals and (1902) 177-227. 


for the rest of the population. In the article, the Cambridge men were 
taken to be representative of the general population. 


Frequency 


60 80 
Height in inches 


Figure 1.2 The heights of 1000 Cambridge men 


Clearly, the geometric distribution in Figure 1.1(a) is not a suitable model 
for the variation in the heights of Cambridge men in 1902, a sample of 
whose heights is represented in Figure 1.2. First, the shapes of the two 
diagrams are different, We would not expect the frequency diagram for 

a sample drawn from a population with a distribution like that in 

Figure 1.1(a) to look like that in Figure 1.2, And secondly, the geometric 
distribution is a discrete distribution — it is a model for a variable which 
can take only the separate, distinct values 1,2,3,...; whereas a height may 
take any value in an interval of values — in this case, between about 

60 inches and 80 inches ~ that is, it is a continuous variable. 


In Chapter D1, the geometric distribution was obtained as a model for the 
variation in the number of rolls of a die needed to obtain a six by assuming 
that each score on a die is equally likely to occur on each roll. There is no 
comparable assumption about heights of men that would allow us to use 
probability theory to develop a model for the variation in the heights of 
men. So a different approach is needed in this case. 


Instead of using theory to develop a model (as we did when considering the 
number of rolls of a die needed to obtain a six), we shall choose a model to 
‘fit’ the data; this model will be one of a collection of standard models. 

If the model that we choose produces a ‘good fit’, then simulations using 
the model should give rise to frequency diagrams which are similar in 
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shape to the frequency diagram for the original data. In particular, the 
two diagrams should have peaks in similar positions and have a similar 
spread of values. So an essential feature of a model for the heights in 
Figure 1.2 is that its probability diagram should be broadly similar in 
shape to the frequency diagram. In addition, the mean value predicted by 
the model should be roughly equal to the mean of the data, and the range 
of values which, according to the model, are likely to occur should 
correspond roughly to the range of values in the data. Finally, the model 
must allow all values in an interval to occur, not just a discrete set of 
values — that is, a continuous model is required. 


Activity 1.1 Properties of the data 


Describe the shape of the frequency diagram in Figure 1.2. Use the 
diagram to estimate very roughly the mean height of the men. Over what 
range are the heights of the men spread? 


Comment 


The frequency diagram is roughly symmetrical about a single peak or 
mode; it is unimodal. The frequencies are low on the left side of the 
diagram, rise steadily to reach a maximum for heights of 69 inches, 

then decrease towards the right. The diagram could be described as 
approximately ‘bell-shaped’. Since the diagram is roughly symmetrical, 
the mean height is approximately 69 inches, the ‘centre’ height. The 
heights of the men range from approximately 61.5 inches to approximately 
77.5 inches. 


So what would the probability diagram of a model for the heights of 
Cambridge men in 1902 look like? Ideally, it should possess all the 
properties of the frequency diagram just described. This suggests that it 
should be ‘bell-shaped’ with a peak at 69 inches and it should allow heights 
between, say, 60 and 78 inches. One possibility is shown in Figure 1.3. 


: > 
60 69 78 
Height in inches 


Figure 1.3 A possible model 


Notice that in this figure the model has been represented by a curve, 
rather than by a series of bars (as in the discrete probability diagram for 

a geometric distribution). It is usual to use a curve when modelling a 
continuous variable. But how is such a model to be interpreted? And how 
can it be used to estimate the proportion of (say) all Cambridge men who 
were over six feet tall, for example, or the proportion who were between 69 
and 71 inches tall? 
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Essentially, areas under the curve are used to represent proportions or 
probabilities. You will see how this is done later in this section when we 
return to the question of the interpretation of the model. The scale on the 
vertical axis, which is omitted in Figure 1.3, is chosen so that areas 
represent proportions. In practice, as you will see later, you do not need to 
worry about including the vertical scale when using a curve as a model: 
often, we shall not even include the vertical axis. The essential point at 
this stage is that the curve should be similar in shape to the frequency 
diagram for the data. 


Activity 1.2 Models for variation 


Frequency 
1s Frequency 4 
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(a) Times between major earthquakes (b) Wing lengths of meadow pipits 


Frequency Frequency 4 
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(c) Durations of eruptions of Old Faithful geyser (d) Weekly earnings of women employees 


Figure 1.4 Frequency diagrams for four data sets 


Figure 1.4 contains frequency diagrams for four data sets. Sources: (a) $237 The Earth: 


‘ . Prater structure, composition and 
Figure 1.4(a) represents the times between ‘major’ earthquakes from oliition, Block 2; 


December 1902 until March 1977. Any earthquake whose magnitude was (b) T. W. Dougall, 


at least 7.5 on the Richter scale or in which over a thousand people were ‘Post-juvenile moult in the 
killed has been included. There were 63 major earthquakes according to Meadow Pipit’, Ringing and 
these criteria, so there are 62 times between earthquakes in the data set. Migration 14 (1993) 137-142. 


Figure 1.4(b) represents the wing lengths of 252 first-year meadow pipits 
netted in southern Scotland in the autumn of 1991. We shall investigate 
these data further in Chapter D4. 


Figure 1.4(c) represents the durations of 106 eruptions of the ‘Old Faithful’ 
geyser in Yellowstone National Park in August 1978. 


Figure 1.4(d) represents the gross weekly earnings (including overtime 
payments) in April 1994 of 150 women working full-time. 
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For each of the four data sets, briefly describe the shape of its frequency 
diagram and sketch a curve which you think might reasonably model the 
variability in the data. Mark on your sketch any key values (as was done in 
Figure 1.3). (Do not worry about the scale on the vertical axis.) 


Comment 
A solution is given on page 45. 


You will probably have noticed that three of these four frequency diagrams 
have one tail of values much more extended than the other. (In fact, the 
frequency diagram for the earthquakes data has no left tail at all.) A data 
set or diagram which is asymmetric because of a long tail is said to be 
skewed. If it has a long right tail, then it is said to be right-skewed; if it 
has a long left tail, it is said to be left-skewed. The terms ‘right-skewed’ 
and ‘left-skewed’ are useful ones for describing the shape of a diagram or 
data set. Figure 1.4(a) and Figure 1.4(d) are right-skewed, and 

Figure 1.4(b) is left-skewed. 


Activity 1.3 Similar models 


Frequency Frequency 4 
30-4 
1000 
20-7 
500 
10>) 
0 
32 34 36 38 40 42 44 46 48 50 0 20 30 40 50 
Chest measurement (inches) Weight (grams) 
(a) Scottish soldiers’ chest measurements (b) Weights of Irish dipper nestlings at 6-8 days 
Frequency Frequency 4 
60>] 
40-7 
2074 
Ly 
—80 -—60 -40 -20 0 20 «40 26 
Radial velocity (kms?) Pengch (runt) 
(c) Radial velocities of stars (d) Lengths of cuckoo eggs 


Figure 1.5 Frequency diagrams for four more data sets 
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Figure 1.5 contains frequency diagrams for four data sets. 
Figure 1.5(a) represents the chest measurements of 5732 Scottish soldiers. 


Figure 1.5(b) represents the weights of 198 Irish dipper nestlings aged 
6-8 days. 


Figure 1.5(c) represents the radial velocities of 80 bright stars in a certain 
region of the sky. (The radial velocity of a star is the velocity with which 
the star appears to be approaching the Earth or, more precisely, the 
component of its velocity towards the Earth. If the radial velocity of a star 
is negative, then the star is moving away from the Earth.) 


Figure 1.5(d) represents the lengths of 243 cuckoo eggs. 


(i) These frequency diagrams are all similar in shape. How would you 
describe their common shape? 


(ii) For each data set, sketch a curve which you think might reasonably 
model the variability in the data, Mark any key values on your 
sketches. (Do not worry about the vertical scale.) 


Comment 


(i) Each of the diagrams is unimodal and roughly symmetrical about its 
peak. In each diagram, the frequencies are low on the left-hand side, 
rise steadily to reach a maximum in the centre, and then decrease 
towards the right. Apart from some jaggedness in Figure 1.5(b) and 
Figure 1.5(d), the diagrams are similar in shape to the frequency 
diagram of heights of Cambridge men in Figure 1.2; they are all 
roughly ‘bell-shaped’. For each data set, the peak is in a different 
position and the spread of the values is different, but all the diagrams 
have the same underlying shape. 

(ii) Sketches of possible models are shown in Figure 1.6. The jaggedness 

in some of the frequency diagrams could simply be due to chance 

variation, In Section 3, you will be asked to investigate whether this 
might be the case, using the statistics software. 


32 40 48 8 28 48 
(a) (b) 
ae 
80 20 40 18 2 26 


(c) 
Figure 1.6 Possible models 


Sources: (a) S. M. Stigler, 
The History of Uncertainty 
before 1900 (Belknap Press, 
Cambridge, MA, 1986); 

(b) John O'Halloran, Patrick 
Smiddy and Barry O’Mahony, 
‘Biometrics, growth and sex 
ratios amongst Irish Dippers’, 
Ringing and Migration 13 
(1992) 152-161; 

(c) R. J. Trumpler and 

H. F. Weaver, Statistical 
Astronomy (University of 
California Press, Berkeley, 
CA, 1953); 

(d) O. H. Latter, ‘The egg of 
Cuculus canorus’, Biometrika 
1 (1902) 164-176. 
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As used here, the word 
‘normal’ is a technical term 
and does not have its 
everyday meaning. 


You do not need to remember 
this formula, 


jeand o are the Greek 
lower-case letters mu and 
sigma. 
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The normal distribution 


The frequency diagrams in Figures 1.2 and 1.5 are all essentially the same 
shape — in each case, a ‘bell-shaped’ curve would seem to be a good model 
for the variation in the data. A model which has this shape is the normal 
distribution. In fact, there is a whole family of normal distributions, all 
with the same basic shape. It is possible to choose a member of this family 
so that its peak is in whatever position is required and with whatever 
spread of values is needed. The equation of a typical normal curve is 

y = f(x), where 


uf 1/xr-p\? 
f(x) - 0 3( = y) (-0o < z < 00). (1.1) 

Notice that (1.1) contains two parameters, 4 and o. Each pair of values for 

and o gives a different member of the family of normal distributions. 

In some texts, symbols other than jz and o are used for the parameters of a 

normal distribution (a and b, for example). However, j: and o are the 

symbols most commonly used, and we shall use them in this course. 


The function f is known as the probability density function of the 
distribution. Note that any value at all for z produces a corresponding 
value for f(x), so the normal distribution is a continuous model. A sketch 
of this curve is shown in Figure 1,7. 


f(z) 


av 


Figure 1.7 A typical normal curve 


The normal curve is symmetrical about the line « = jz, so the parameter jz 
gives the location of the ‘centre’ of the curve. Since the heights of 
Cambridge men in Figure 1.2 are distributed roughly symmetrically about 
69 inches, we might use a normal distribution with parameter j. = 69 to 
model the variation in the heights of Cambridge men. Similarly, a normal 
distribution with 4. = 40 might be used to model the variation in Scottish 
soldiers’ chest measurements shown in Figure 1.5(a). 


SECTION 1 MODELS FOR VARIATION 


The parameter o governs the spread of the values that, according to the 
model, are most likely to occur: in general, for fixed horizontal and vertical 
scales, a small value of o produces a tall narrow ‘bell’, and a large value of 
o produces a short wide ‘bell’, Two examples are shown in Figure 1.8. The 
roles of the two parameters and the reason for the use of the symbols 4 
and o are discussed further in Section 2. 


RY 


so 100 120 40 60 80 100 120 140 160 # 
(a) Smaller o (b) Larger o 


Figure 1.8 Two normal curves (same vertical and horizontal scales) 


In Section 2, we shall discuss how to choose values for jz and o in order to 
select an appropriate normal distribution as a model for the variation 
observed in a data set; and in Section 3, you will have the opportunity to 
fit normal models to a number of sets of data using the statistics software. 


A normal curve is symmetrical about its peak and falls away on either 
side, gradually approaching the z-axis. Note, however, that the curve 
never actually meets the x-axis, unlike the sketches of possible models 
drawn in Activity 1.3. In this model, any number between —oo and +00 is 
possible, although values far from the peak are very unlikely. You might 
wonder whether it is reasonable to model measurements such as height, 
length and chest size, which cannot be negative, using a distribution which 
allows negative values. The point to remember is that a normal 
distribution is only a model; in practice, when fitting a normal model to 
such data, the proportion of values that, according to the model, are 
negative is so small that it may be ignored. You can see this for yourself in 
the computer session in Section 3, 
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Using a curve to model variability 


To see how a normal curve can be used to model variability, we return to 
the data on the heights of 1000 Cambridge men which were represented in 
the frequency diagram in Figure 1.2. This diagram is repeated below as 
Figure 1.9(a). 


Frequency 4 
200-4 
160+ 
1204 


804 


% 80 

(a) Height in inches 
Proportion 4 
0.20-4 

0.167} 

0,124 

0.087) 


0,04>) 


Li 


60 65 70 75 


(b) Herein toons 


Figure 1.9 (a) A frequency diagram (b) A histogram 


Figure 1.9(b) shows a histogram for the same data. While in a frequency 
diagram the heights of the bars represent the various frequencies, in a 
histogram it is the areas of the bars which are proportional to the 
frequencies. To obtain the histogram in Figure 1.9(b) from the frequency 
diagram in Figure 1.9(a), the scale on the vertical axis was adjusted so 
that the sum of the areas of all the bars is equal to 1. In all other respects, 
the two diagrams are identical. Since the total area of the bars in the 
histogram is equal to 1, this means that, for example, the area of the bar 
centred on 70 inches, which represents heights between 69.5 and 70.5 
inches, is equal to the proportion of the 1000 men who were between 
69.5 and 70.5 inches tall. In Figure 1.9(a), the height of this bar is 139, 
since 139 of the 1000 men had heights in this range; in Figure 1.9(b), the 
area of the corresponding bar is 0.139, since this is the proportion of men 
whose heights were in this range. 


Suppose now that an appropriately chosen normal curve is superimposed 
on the histogram in Figure 1.9(b); this is shown in Figure 1.10. (You will 
learn how to choose an appropriate normal curve in the next section.) 


SECTION ! MODELS FOR VARIATION 


60 65 70 80 
Height in inches, 


Figure 1.10 The histogram and a fitted normal curve 


The total area between the normal curve and the z-axis may be obtained 
as follows. 


The area between the normal curve and the z-axis from x = —N to «= N 
is given by 
N 
| f@rae, 
-N 


where f(x) is given by formula (1.1). This integral can be calculated for 
any N > 0, and its value approaches a limiting value as NV becomes larger 
and larger. The total area between the normal curve and the z-axis is 
given by this limiting value. We write 


area = [ f(e)de = Jim ii flo} ae. 
co 00 JN 


For any normal curve this area is always equal to 1, no matter what are Some information about 

the values of the parameters jz and o. So the total area under any normal _ several results involving 
curve is equal to the total area of the bars in a histogram, provided the integrals with infinite limits 
histogram is drawn so that the total area of the bars is equal to 1. that are stated without proof 


in this chapter are included in 
As you can see, the area under the part of the curve in Figure 1.10 from an appendix (page 43). 
69.5 to 70.5 is approximately equal to the area of the corresponding bar in 
the histogram. These areas are shown in Figure 1.11. 


69.5 70.5 


Figure 1.11 The area under the curve and the area of the bar 
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The shaded area under the curve in Figure 1.11 is, in fact, equal to 0.141 
to three decimal places. (In Section 3, you will learn how to use the 
statistics software to find this area.) So the area under the curve between 
69.5 and 70.5 is approximately equal to the proportion of men in the 
sample who were between 69.5 and 70.5 inches tall — 0.141 in the model 
compared with 0.139 in the sample. 


The histogram represents the heights of just 1000 Cambridge men in 1902. 
If we assume that these were a sample from the population of all men in 
Cambridge, then the normal curve provides a model for the variation in 
the heights of all Cambridge men in 1902. The area under the curve 
between the two heights 69.5 inches and 70.5 inches represents the 
proportion of all Cambridge men that were, according to the model, 
between 69.5 and 70.5 inches tall. 


We can generalise the interpretation of areas under the normal curve as 
follows. An area under the curve between two heights, a inches and 

b inches say, represents the proportion of all Cambridge men that were, 
according to the model, between a inches and b inches tall. This is 
illustrated in Figure 1.12. 


Proportion of all 
Cambridge men 
between 
a and b inches 
tall 


Figure 1.12 The model 


The shaded area could be found by evaluating the integral 


b 
f(v) de, 
where f(x) is given by formula (1.1). 


Since the area represents the proportion of Cambridge men that were 
between a and 6 inches tall, it follows that if a man had been selected at 
random from all the men in Cambridge in 1902, then the probability that 
his height would have been somewhere between a inches and 6 inches is 
also given by this area. This gives us a second interpretation of the area. 
So, for instance, if a man had been selected at random from all Cambridge 
men, then the probability that his height would have been between 

69.5 and 70.5 inches is given by the area under the curve between 69.5 
and 70.5. 
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Activity 1.4 Interpreting the model 


Figure 1.13 contains three sketches of the particular normal curve used to 
model the heights of all Cambridge men in 1902. For each sketch, describe 
in words what the shaded area represents. 


69 71 


. Sowa 


(a) 


Figure 1.13 Three areas under a normal curve used to model the heights of Cambridge men in 1902 


Comment 
A solution is given on page 45. 


You might wonder why we should use a model to estimate the proportion 
of Cambridge men between, say, 69 and 71 inches tall. Surely it would 
have been better to measure all the men and find the actual proportion 
whose height was in this range? In theory, it would have been better, but 
in practice, this would not have been a realistic course of action to take. 
The work involved in measuring the heights of 1000 men must have been 
very considerable, To measure the height of every adult man in Cambridge 
would have been a yery difficult, if not impossible, task. Consider the 
difficulties involved in such an undertaking! How could you be sure that 
no one had been missed? 


This sort of situation is a common one. Information about an entire 
population is needed (perhaps for developing marketing or manufacturing 
plans, or for making decisions about the future provision of health care or 
education, or perhaps simply to gain knowledge), but it is impractical or 
too costly to examine (or interview or measure, as appropriate) every 
member of the population. A practical approach is to examine a sample 
drawn from the population. 


If the sample is suitably selected (ideally, it should be a random sample 
that is, a sample drawn in such a way that every member of the population 
has an equal chance of being chosen), then we can use statistics to infer 
information about the population from what we find out about the sample. 
And in order to make inferences about the population, we frequently need 
to develop a model for the variability in the population as a whole. 


What we have been doing is seeking a model for the variation in the 
heights of all Cambridge men which fits the data available for the sample 
of 1000 men. This model can then be used to estimate the proportion of 
all Cambridge men whose height is in any given range. Such proportions 
may be found by calculating appropriate areas under the curve used to 
model the variation in heights. 
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A little history 


The normal distribution is one of the most frequently used models for 
variation; its specification given in (1.1) was obtained first by Abraham 
de Moivre (1667-1754). He derived it as an approximation while 
endeavouring to compute probabilities in various games of chance. 
After the existence of a belt of asteroids between Mars and Jupiter was 
discovered at the beginning of the 19th century, the study of the 
motions of the asteroids was the subject of considerable scientific 
activity. One of the very early applications of the normal distribution 
was in describing the variability in errors of measurement of the 
motions of such asteroids. 


The normal distribution is also known as the Gaussian distribution, 
after the German mathematician and astronomer Karl Gauss 
(1777-1855), who played a leading role in demonstrating the usefulness 
of the distribution as a model for measurement errors in general in the 
physical sciences. Since then, the distribution has been used 
successfully to model variability in many phenomena in other fields of 
study. Indeed, for a time in the 19th century, there was a school of 
thought which maintained that the normal distribution was the only 
model necessary to describe variability, and a great deal of effort was 
put into trying to prove that this was the case, 


Although the normal distribution is a very useful model, it is not 
universally applicable. You have already seen this; in Activity 1.2, you 
met four examples where a bell-shaped curve would not have been an 
appropriate model. But as you will see in Chapter D3, in certain 
circumstances the normal distribution may still be used in a statistical 
investigation, even when it is not a suitable model for the variability in the 
available data, and this is the reason for its importance in statistics. 


Activity 1.5 Making notes, preparing a presentation 


The idea of using a curve to model variation may be new to you. Spend a 
few minutes thinking about what you now understand about how a normal 
distribution can be used to model the variation observed in data, and write 
down your ideas in your own words. 


Towards the end of this chapter you will be invited to consider how you 
would prepare and give a presentation on modelling variation to other 
students (see Activity 7.1). As you continue to work through the sections, 
think about how you could use your own notes to help you to prepare such 
a presentation, For example, what do you consider to be the key point(s) 
from this section that you want to record and so include as part of your 
presentation? Are there any diagrams you intend to include? You may 
find it helpful to record your ideas on Learning File Sheet 2. 


Preparing a presentation encourages you both to be self-evaluating and to 
give and receive constructive feedback; presenting work encourages the 
need for clear and logical structure, as well as effective use of diagrams and 
other visual aids, 
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Summary of Section 1 


In this section, you have seen how a curve may be used to model the 
variation in various quantities, such as men’s heights, women’s earnings, 
the times between major earthquakes and the wing lengths of meadow 
pipits. The frequency diagrams for several of the data sets had the same 
‘bell-shape’. A distribution was introduced as a model for variation of this 
sort — the normal distribution. You saw how areas under a normal curve 
can be used to answer questions about the population being modelled. 


The equation of a normal curve contains two parameters, 1 and o. 

In order to select a particular model, values for these parameters must be 
chosen. In the next section, we discuss what these parameters represent 
and how suitable values for them can be estimated from the available data. 
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2 Choosing a normal model 


# (mu) is the Greek letter ‘m’ 
and o (sigma) is the Greek 
letter ‘s’. There is therefore a 
first-letter link (albeit in a 
different alphabet) between 
the symbol for the parameter 
and what it stands for. 


In general, letters from the 
Roman alphabet are 
commonly used for sample 
Statistics, and letters from the 
Greek alphabet are used for 
population parameters. 
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In Section 1, you saw that the equation of a normal curve contains two 
parameters, denoted 1 and o. In order to choose a specific normal curve to 
model the variability in some quantity — such as, for example, the heights 
of Cambridge men or the weights of Irish dipper nestlings — you need to 
specify values of these parameters. The purpose of this short section is to 
discuss what they represent and how to select values for them. 


In Section 1, it was observed that the parameter j: provides the location of 
the ‘centre’ of the curve, and the parameter o governs the spread of values 
that are most likely to occur. In fact, the parameter jis the mean of the 
distribution and o is a measure of spread called the standard deviation. 
The symbols 1 and o are used generally to denote the mean and standard 
deviation of a distribution. In practice, probability distributions are used 
to model the variation in populations, and so the symbols jand @ are also 
used to denote the mean and standard deviation of a population. 


The use of the symbols j and o for the mean and standard deviation of 
both distributions and populations does not usually lead to any confusion, 
because of the close connection between a population and a probability 
distribution used to model the variation in the population. In fact, the 
mean of a probability distribution is sometimes called the population mean 
and the standard deviation of a probability distribution is sometimes 
called the population standard deviation. The population mean wand the 
population standard deviation o are called population parameters. 


One of the data sets in Section 1 consisted of the heights of 1000 
Cambridge men in 1902. Suppose that we wish to use a normal curve to 
model the variability in the heights of all Cambridge men in 1902. If the 
model is to reflect accurately the distribution of heights in the population, 
then the parameter 1 should be equal to the mean height of all Cambridge 
men ~ the population mean ~ and the parameter o should be equal to the 
standard deviation of the heights — the population standard deviation, 

So we ought to use the population mean and the population standard 
deviation as values for the parameters of the normal model, 


But we do not know the mean and standard deviation of the population. 
However, we can use the sample of 1000 heights drawn from the 
population to calculate the corresponding sample statisti the sample 
mean, which is denoted %, and the sample standard deviation, denoted s. 
The sample mean F may then be used as an estimate for the population 
mean, and the sample standard deviation s as an estimate for the 
population standard deviation. The approach that we shall adopt in 
Section 3, when we use the statistics software to fit a normal curve to data, 
is to use the values of the sample mean Z and the sample standard 
deviation s to estimate the population mean jz and the population 
standard deviation o, respectively. So the values of the parameters of the 
normal model that we fit will be given by 7 and s. 


In this section, we shall discuss briefly how the population parameters 4: 
and o are defined, and how the sample statistics F and s are calculated. 
The population mean jz and the sample mean Z are discussed in 
Subsection 2.1; the population standard deviation o and the sample 
standard deviation s are discussed in Subsection 2.2. This section contains 
quite a lot of detail that you are not expected to remember. For instance, 
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in this course, you will not be expected to calculate the mean and standard 
deviation of a probability distribution. When these are required 

(in Chapter D3, for example), you will be given their values. But after 
working through this section, you should make sure that you understand 
the distinction between population parameters such as yz and o and sample 
statistics such as Z and s, and that you can calculate the mean and 
standard deviation of a small data set using a calculator. In Section 3, you 
will learn how to use the statistics software to calculate the sample mean 
and sample standard deviation. 


2.1 The mean 


Many of the ideas in this subsection will be familiar to you. For instance, 
you should already be familiar with the calculations involved in finding 
the mean of a sample of data (perhaps from studying MU120 or the 
preparatory materials). And the mean of a discrete distribution (such as 
the geometric distribution) was introduced in Section 4 of Chapter D1. 
Before discussing how the mean of a continuous distribution such as the 
normal distribution is defined, we shall review briefly how the mean of a 
sample of data is calculated. 


The sample mean 

Suppose that a sample of n observations is taken from a population. 

(For example, these might be the heights of n men.) We shall denote these 
observations by x1, 22,...,@n- The sample mean Z, which is defined to be 
the sum of the observations divided by the number of observations, is 
given by the formula 


1 1 
Fane (Gash Prt inerhahens Dee 


Sometimes, many of the recorded values occur several times. In this case, 
the data are often placed in a table showing the number of times that each 
different value was observed ~ its frequency. For example, Table 2.1 
contains the data on the heights of 1000 Cambridge men which were 
represented in Figure 1.2. 


Table 2.1 Heights of 1000 Cambridge men 


Height in inches Frequency 


62 3 
63 20 
64 24 
65 40 
66 85 
67 122 
68 139 
69 179 
70 139 
7 107 
72 55 
73 47 
74 22 
75 12 
76 5 
77 1 


You might like to think about 
why the two formulas are 
equivalent. 
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To find the mean height when data are given in this form, we use the 
equivalent formula 


Le 
Ee 
ny 


i= Ahi + aah +--+ 2h) 


where a ,22,...,2, are the different values observed in the sample and 
Si, fo, ---, fy ave their respective frequencies; the sample size n is equal to 


fit fo+-::+ fx, the sum of the frequencies. 

For the sample of 1000 heights of Cambridge men, 
1 
1000 
So the sample mean is approximately 68.9 inches. The term in brackets is 
the sum of the heights of all 1000 men: 3 men were 62 inches tall, so 62 is 
included 3 times (62 x 3); 20 men were 63 inches tall, so 63 is included 

20 times (63 x 20); and so on. 


(62 x 3-+63 x 20+++-+77 x 1) = 68.872. 


z= 


In fact, the heights given in Table 2.1 were recorded rounded to the 
nearest inch, so that a man recorded as 62 inches tall may have been as 
much as half an inch shorter or taller than 62 inches. All men who were at. 
least 61.5 inches tall but less than 62.5 inches tall were grouped together 
and recorded as being 62 inches tall. Similarly, all men between 68.5 and 
69.5 inches tall were grouped together and recorded as 69 inches tall. So it 
is possible, even likely, that the value of the sample mean calculated above 
is not exactly equal to the actual mean height of the 1000 men. However, 
given the way the data have been recorded, it is the best that we can do: 
the midpoint of each interval of heights has been used instead of the exact 
heights of the men in each group. 


Data are often summarised and reported by grouping the values and 
recording the frequencies for each group. In this case, the midpoint of each 
interval is used in calculating statistics such as the sample mean. This was 
done implicitly when calculating the mean height of the Cambridge men. 


Activity 2.1 Irish dipper nestlings 


Table 2.2 Weights of Irish dipper nestlings 


Weight in grams Frequency 


9-11 1 
11-13 8 
13-15 3 
15-17 6 
17-19 5 
19-21 12 

- 21-23 13 
23-25 18 
25-27 27 
27-29 18 
29-31 22 
31-33 20 
33-35 3 
35-37 17 
37-39 8 
39-41 7 
41-43 4 
43-45 4 
45-47 2 
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The data which were used to draw the frequency diagram in Figure 1.5(b) 
are given in Table 2.2. These are the weights of 198 Irish dipper nestlings 
aged 6-8 days. What values would you use to calculate the mean weight of 
these nestlings? (Do not bother to go on to calculate the mean: leave this 
sort of lengthy calculation for the computer to do. You will learn how to 
do this in Section 3.) 


Comment 
The midpoints of the intervals are 10, 12, 14,...,46, so these are the values 
that are to be used to calculate the sample mean. The sample mean is 
given by 

ul 


7 S198 


(10 x 1412 x 8+++- +46 x 2). 


The population mean 


In Section 4 of Chapter D1, you were introduced to the idea of the mean 
of a probability distribution or the mean of a random variable X. In the 
situation discussed there, X was the number of rolls of a die needed to 
obtain a six, and you saw that this random variable has a geometric 
distribution. You also saw that, for a discrete distribution such as the 
geometric distribution, the mean of the distribution, which is denoted by 
the Greek lower-case letter j, is given by the formula 


Dix P(x =3), (2.1) 
7 
where the summation is over all values j that X can take. The probability 
distribution provides a model for the variability in the number of rolls of 
the die needed to obtain a six; so the mean of the distribution is the mean 
number of rolls predicted by this model. 


Since a probability distribution may be used to model the variability in a 
population, the mean of the distribution is also sometimes called the 
population mean. (The population in the example above consists of the 
possible numbers of rolls of a die needed to obtain a six.) 


For a discrete model, the population mean is given by formula (2.1). But 
clearly this is not appropriate for a continuous model, such as the normal 
distribution, in which all real values in an interval are possible. For a 
continuous model, variability is described by a curve, instead of a discrete 
set of probabilities, and areas under the curve represent probabilities. 


The formula for the mean of a continuous distribution, which corresponds 
to formula (2.1) for the mean of a discrete distribution, is 


mean =/ xf(x) da, (2.2) 
where f is the probability density function of the distribution (that is, 

y = f(z) is the equation of the curve used as a model). Formula (2.2) is 
the continuous analogue of formula (2.1): an integral replaces the sum, and 
f(x) replaces the discrete probabilities P(X = j). 


See Chapter D1, 
formula (4.4). 


See the appendix to this 
chapter for further details. 


The word ‘mean’ is used twice 
in this phrase in two slightly 
different roles: deviations 
from the mean (seen as a 
place), and then the mean of 
these deviations, to be 
calculated. 
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For a normal distribution with parameters j: and o, the integral on the 
right-hand side of formula (2.2) evaluates to the parameter yz contained in 
the equation of the normal curve. So the mean of the distribution is equal 
to the parameter jz. (You might have expected this result: a normal curve 
is symmetrical about the value x = ju, so the mean is 1.) So the use of the 
symbol ju for the parameter of a normal distribution reflects the role of the 
parameter in the model. 


In general, given a sample of data, the sample mean Z is used to estimate 
Ht, the mean of the population from which the sample was drawn. So if 
we have a sample from a population and we wish to choose a normal 
distribution to model the variation in the population, then we shall use 
the sample mean Z for the parameter ju in the model. For the heights of 
Cambridge men, the sample mean is approximately 68.9 inches, so we 
might well choose a normal distribution with parameter j: = 68.9 to model 
the variation in the heights of all Cambridge men. 


2.2 The standard deviation 


There are several measures of spread which may be used to summarise the 
variability in a sample of data. If you have studied MU120, then you will 
already be familiar with the range, the interquartile range, the relative 
spread and the standard deviation as different possible measures of spread. 
As indicated at the beginning of this section, the parameter o of a normal 
distribution is the standard deviation of the distribution, so in this 
subsection we shall restrict our attention to this particular measure of 
spread. We shall discuss how it is defined and calculated, both for a 
distribution and for a sample of data. 


The population standard deviation 


Consider first a population with mean jz. The population standard 
deviation is a measure of how widely the values in the population are 
spread about the mean ji; it is defined in terms of the deviations of values 
from ju. The deviation of a value x from the mean pu is « — 1; so the 
squared deviation is written (« — )*. The population standard deviation 
is defined most easily in terms of a related measure called the population 
variance. The population variance is defined to be the mean squared 
deviation from the mean; that is, it is the average (i.e. the mean) value of 
(x — )® for the whole population. 


You have just seen that, for a continuous model, the average or mean of 
the values x in a population is found by integrating the product xf (x) over 
all possible values x, where f is the probability density function of the 
distribution which describes the variation in the population. The mean of 
the squared deviations (x — j.)? is found in a similar way, but this time we 
integrate the product (2 — 1)? f(x) over all possible values of x: 


es 
population variance = " (@ — yu)? f(x) de. (2.3) 
oo 
The population standard deviation is defined to be the square root of 
the population variance. 


For a normal distribution with parameters 4: and o, with probability 

density function given by formula (1.1), the integral on the right-hand side 

of formula (2.3) evaluates simply to 07. So, since the standard deviation is See the appendix. 
the square root of the variance, the standard deviation of the distribution 

is equal to a, as you were told earlier. As already mentioned, the symbol a 

is used generally to denote the standard deviation of a population or of a 

probability distribution. So the use of the symbol o for the second 

parameter of a normal distribution also reflects the role of the parameter 

in the model. 


The sample standard deviation 


Now suppose that a sample x,22,...,2,, of size n has been obtained from 
a population with unknown mean j and unknown standard deviation oc. 
You have seen that the sample mean Z can be used to estimate the 
population mean ji. Correspondingly, the sample standard deviation, which 
is denoted by s, can be used to estimate the population standard 
deviation. But how is the sample standard deviation defined? 


Since the population variance is defined to be the mean squared deviation 
from the population mean 1, by analogy we could define the sample 
variance to be the mean squared deviation from the sample mean 7. 
Replacing y by % and dividing the sum of the squared deviations by n to 
find the mean squared deviation, this would suggest the formula 


, yp ae 
variance = — vj- Tt). 
n >a ) 


Hence, taking the square root, 


standard deviation = 


This is the formula for the standard deviation that was introduced in 
MU120. However, in practice, this formula tends to underestimate the 
population standard deviation o. If we are fitting a normal model to some 
data, then we would not want to choose a model which underestimates the 
variability in the population from which the data were drawn. So we would 
not want to use precisely this formula to choose the value of the parameter 
oa. It can be shown, though we shall not try to justify it, that this tendency 
to underestimate the population standard deviation a can be corrected by 


using the divisor n — 1 rather than n. (When n is large, the difference ‘The proof of this result 

between the values obtained using a divisor n and a divisor n — 1 is requires ideas and techniques 

negligible.) This results in the following definitions. which are beyond the scope of 
this course. 


The sample variance, denoted by s”, is defined to be 


Seek 
28 ac eseNl 
s Sata 2)*; 


and the sample standard deviation, denoted by s, is given by 


n 


(2.4) 
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The formula for the sample standard deviation looks complicated, but the 
idea behind it is straightforward and, in practice, your calculator or your 
computer will do the actual calculations involved. Most calculators provide 
two different versions of the standard deviation, One popular brand has 
two options labelled o,,_; and o,,; another has two options labelled Sx and 
ox. These options correspond to the two different formulas: in each case, 
the first uses the divisor n — 1 and the second the divisor n. If your 
calculator has options for the two different versions of the standard 
deviation, then make sure you know which option is which. Two exercises 
are included at the end of this section so that you can check that you can 
use your calculator to find the sample mean and sample standard deviation 
(that is, with divisor n — 1) of a data set. 


Both formulas for the standard deviation are in common use, but they are 
used in different circumstances. It is usual to use the sample standard 
deviation s with divisor n — 1 as defined above when estimating a 
population standard deviation o. So, when a normal distribution is used to 
model the variation in a population, the sample standard deviation s with 
divisor n — 1 is used for the parameter o in the model. For the heights of 
Cambridge men, the sample standard deviation is approximately 

2.57 inches, so we might well choose a normal distribution with parameter 
ao = 2.57 to model the variation in the heights of all Cambridge men. 


In this course, we shall always use the divisor n — 1 when calculating the 
standard deviation of a sample of data. The statistics software OUStats 
also uses the divisor n — 1. 


Summary of Section 2 


In general, given a sample of data from a population, the sample mean 7 
is used to estimate the population mean jp, and the sample standard 
deviation s is used to estimate the population standard deviation o. 
This is summarised in the table below. 


Sample statistic Population parameter 
bp is used to estimate bw 
8 is used to estimate o 


The sample statistics and s are defined as follows. 


For a sample of n observations 2, 2,...,#,, the sample mean F is given by 


Pee dine 
= a +a +++ +a) = Yay 
t=1 


al 


and the sample standard deviation s is given by 
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More generally, if 21, 22,...,, are the different values observed in a 
sample of n observations, and f;, fo,-.-, fx are their respective frequencies 
(so f, + fo +++: + fe =n), then the sample mean 7 is given by 


Eo 


k 
= dah + 2ofot-++ + 2nfi) = Lath 
n n = 


and the sample standard deviation s is given by 


7 


i=l 


When choosing a normal distribution to model the variation in a 
population, the values of two parameters have to be specified. These 
parameters are the mean and standard deviation of the normal 
distribution. When fitting a normal model to data, the values of the 
sample mean Z and the sample standard deviation s are used for the 
parameters of the model. 


Exercises for Section 2 


It is not always convenient to use a computer to calculate statistics such as 
the sample mean and the sample standard deviation; and it may not be 
worthwhile switching on a computer if the data set is fairly small. So it is 
important for you to be able to perform these calculations on your 
calculator. If your calculator has options for the sample mean and the 
sample standard deviation, then make sure that you can use them to carry 
out the following exercises. If your calculator does not have these 
statistical facilities, then do the calculations by completing a table with 
columns headed «;, 2; — %, (a; — F)*. The sample mean % can be found by 
summing all the numbers in the first column and dividing by n. To find 
the sample standard deviation, add up the numbers in the third column to 
find 33(a; — @)?, divide by n — 1, and finally find the square root. 


Exercise 2.1 How much mail? 


The number of items of mail delivered on a Monday morning to each of 
five houses chosen at random from those in a large estate were as follows. 


Pe cies ae 
(a) Find the sample mean and the sample standard deviation. 
(b) If there are 3000 houses on the estate, what is your estimate of the 
total number of items of mail delivered to the estate on that morning? 
Exercise 2.2 Earnings of mechanical engineers 


The gross weekly earnings (in pounds) in 1995 of a sample of six 
mechanical engineers were as follows. 


310 635 464 520 381 732 


Find the sample mean and the sample standard deviation. 


These are invented data, but 
they are based on statistics 
published in the 1995 New 
Earnings Survey. 


3 Fitting a normal model 
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To study this section, you will need access to your computer, together with 
the statistics software and Computer Book D, and an audio tape player 
together with Audio Tape 3. 


In Section 1, the normal distribution was introduced as a model for the 
variation in various phenomena, including the heights of Cambridge men, 
the chest measurements of Scottish soldiers and the lengths of cuckoo eggs. 
You saw that a particular normal model is chosen by fixing the values of 
two parameters, j: and a. In Section 2, we discussed what these parameters 
represent and how, given a sample of data, their values may be estimated. 
In this section, the software OUStats will be used to fit normal models to a 
number of data sets. You will need your computer throughout this section; 
it has been divided into three subsections to provide suitable break points 
in case you are unable to complete all the activities in one session. 


Before you begin to use OUStats, look back to Learning File Sheet 1 where 
you identified what you might do differently in becoming familiar with a 
software package. 


3.1 Introducing OUStats 


In this subsection, the main facilities of OUStats are introduced, using an 
audio tape session. The data on the heights of Cambridge men are used to 
illustrate the fitting of a normal model. You will also see how the model 
can be used to answer questions such as ‘what proportion of all Cambridge 
men in 1902 were between 69.5 and 70.5 inches tall?’. 


Refer to Computer Book D for the work in this subsection. 


You will need to listen to Audio Tape 3, Band 3, ‘Introducing OUStats 
for MST121’, during the computer work. 


3.2 Is a normal model a good fit? 


In this subsection, you will be asked to explore several data sets and, in 
each case, to investigate whether a normal distribution provides a good 
model for the variation observed in the data. 


Refer to Computer Book D for the work in this subsection. 


3.3 Printing with OUStats 


At some stage in your study of this block, you will need to print output 
from OUStats. For instance, you may need to include diagrams produced 
by OUStats in your answers to TMA questions. This subsection contains 
instructions for printing output from OUStats. It would be a good idea to 
follow these instructions in order to check that you can use your printer to 
obtain such output. 


@) Refer to Computer Book D for the work in this subsection. 


SECTION 3 FITTING A NORMAL MODEL 


Summary of Section 3 


In this section, you have begun to use OUStats, the data analysis part of 
the statistics software. You have used OUStats in the following ways. 


oo 9 099090 


To obtain summary statistics for data. 

To obtain frequency diagrams for data. 

To fit a normal curve to data. 

To calculate areas under a normal curve. 

To generate random samples from a normal distribution. 


To print output. 


You have also used your experience from one software package (Mathcad) 
to help in learning and using a different one (OUStats). : 
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To study this section, you will need access to your computer, together with 
Computer Book D and the statistics software. 


In this section, you will have the opportunity to investigate a large data 
set on the heights of fathers and their sons. The data were collected in the 
1890s as part of a larger study into the inheritance of characteristics in 
humans. The data are of historical interest and have been the subject of 
much study since they were published. 


Refer to Computer Book D for the work in this section. 


Summary of Section 4 


In this section, you have investigated a large data set on the heights of 
fathers and sons. You have used many of the facilities of OUStats that 
were introduced in Section 3. In addition, you have used OUStats to 
obtain a scatterplot for paired data. 


To study this section, you will need access to your computer, together with 
Computer Book D and the statistics software. 


In Section 1, the normal distribution was introduced and some basic 
properties of normal distributions were described: for example, the position 
of the ‘centre’ of a normal curve is given by the parameter j1, which is the 
mean of the distribution, and the spread of values is governed by the 
parameter ¢, the standard deviation of the distribution. In Sections 3 

and 4, you used the statistics software to fit normal models to data and to 
do calculations for normal distributions. In this section, you are invited to 
use the software to explore further the properties of normal distributions. 


Refer to Computer Book D for the work in this section. 


Summary of Section 5 


Several general properties of normal distributions have been illustrated in 
the activities in this section. These and other properties of normal 
distributions are summarised in the final section of this chapter. 


C8) 
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The normal distribution underpins much of the work that you will be 
doing in Chapters D3 and D4. The purpose of this section is to bring 
together the main properties of normal distributions that have been 
discussed and explored in this chapter. We begin by summarising the 
properties of normal distributions that were discussed in Sections 1 and 2. 


The normal distribution was introduced as a model for the variation 
observed in a number of samples of data, including the heights of 
Cambridge men, the lengths of cuckoo eggs and the weights of Irish dipper 
nestlings. In fact, there is not just one normal distribution but a whole 
family of closely related normal distributions, each being specified by two 
parameters, most commonly denoted by yz and o. 


The equation of a normal curve is given by 


Shlet a=#) = 
fa) = Faxes 3( = ) (-00 < x < 00). 
A typical normal curve is shown in Figure 6.1. 
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Figure 6.1 A typical normal curve 


A normal curve is symmetric about its peak, which occurs at x = ju; so the 
mean of the distribution is equal to the parameter jz. And the position of 
the peak of a normal curve along the x-axis depends on the value of its 
mean, the parameter j. Normal curves with different means but with the 
same value of o are identical in shape but their peaks are located at. 
different positions on the z-axis. Figure 6.2 shows sketches of the three 
normal curves from Activity 5.1 in Computer Book D; the value of the 
parameter o is 1 in each case. 
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Figure 6.2 Three normal curves each with o = 1 


SECTION 6 PROPERTIES OF NORMAL DISTRIBUTIONS 


The parameter o of a normal distribution is the standard deviation of the 
distribution and is a measure of the spread of values in the distribution. 
Roughly speaking, a small value of o produces a tall narrow ‘bell’ and a 
large value of o produces a short wide ‘bell’. (Although, of course, whether 
a ‘bell’ looks short and wide or tall and narrow will depend on the scales 
used on the axes.) Two normal curves with the same mean but different 
standard deviations are sketched in Figure 6.3; the same vertical scale was 
used to produce the two sketches. 
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(a) o=6 (b) o = 20 
Figure 6.3 Two normal curves (same scales) each with = 100 


You have seen that an area under a normal curve between x = a and x = b 
represents the proportion of the population being modelled which have 
values between a and b, When modelling the heights of Cambridge men in 
1902, for example, the area between x = 69 and x = 71 under the normal 
curve used to model heights represents the proportion of all Cambridge 
men who were between 69 and 71 inches tall. This area is the shaded area 
in Figure 6.4. 


69 71 
Figure 6.4 An area under a normal curve 


An alternative interpretation of this area is as a probability. If a man had 
been chosen at random from all Cambridge men in 1902, then this area 
represents the probability that his height would have been between 69 and 
71 inches. You found a number of areas under normal curves in Section 4 
when investigating the heights of fathers and their sons. 


In Section 5, after finding areas under a small number of normal curves, 
you were invited to make hypotheses about areas for all normal 
distributions, whatever the values of 4 and o. The results of Activities 5.1 
to 5.4 in Computer Book D may be summarised as follows. 
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Suppose that a normal distribution is used to model the variation in a 

population. Then according to the model: 

© approximately 68.3% of the population are within 1 standard 
deviation of the mean (that is, between  — o and p+); 

© approximately 95.4% of the population are within 2 standard 
deviations of the mean (that is, between js — 20 and yz + 2c); 

© almost all the population — about 99.7% — are within 3 standard 
deviations of the mean (that is, between  — 3 and p+ 3a). 

Note that these results hold true whatever the values of the mean j: and 

the standard deviation o of the distribution. The results are illustrated in 

Figure 6.5. 
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Figure 6.5 Areas under normal curves 


Example 6.1 Heights of Cambridge men 


The normal distribution used to model the heights of Cambridge men in 
1902 has mean ps = 68.9 and standard deviation o = 2.57. According to 
this model, within what range were the heights of almost all Cambridge 
men — that is, about 99.7% of them? How does this compare with the 
sample of heights? 


Solution 


Almost all Cambridge men — about 99.7% of them — were between 
68.9 — 3 x 2.57 and 68.9 + 3 x 2.57 inches tall, that is, between 
approximately 61.2 and 76.6 inches. 


So the model predicts that about three in a thousand men will be either 
shorter than 61.2 inches or taller than 76.6 inches. As you saw, for the 
sample of 1000 men from this population, only one man’s height was 
outside the range 61.2 inches to 76.6 inches — no man was less than 

61.5 inches tall, one man was recorded as 77 inches tall and no man was 
taller than this. So ‘almost all’ the men in the sample were within three 
standard deviations of the mean height. The normal model reflects quite 
well the proportion of men in the sample who were unusually short or tall. 


SECTION 6 PROPERTIES OF NORMAL DISTRIBUTIONS 


Activity 6.1 Lengths of cuckoo eggs 


In Activity 3.5 in Computer Book D, a normal distribution with mean 
= 22.4 and standard deviation o = 1.08 was fitted to the data on the 
lengths of cuckoo eggs. According to the model, within what range are the 
lengths of almost all cuckoo eggs (that is, 99.7% of them)? How does this 
compare with the sample of lengths? 


Comment 


A solution is given on page 45, Int THE JURASSIC PERIOD, PTERODACTAS 
LAIDTHEIR. GGS IN CUCKOOS! NESTS. 


Further similar results for a population modelled by a normal distribution 
with mean j and standard deviation o were illustrated in Activities 5.5 
to 5.7 in Computer Book D. These may be summarised as follows: 


© approximately 90% of the population are within 1.64 standard 
deviations of the mean, that is, between p — 1.640 and p+ 1.640; 

© approximately 95% of the population are within 1.96 standard 
deviations of the mean, that is, between  — 1.960 and p + 1.960; 

© approximately 99% of the population are within 2.58 standard 
deviations of the mean, that is, between ps — 2.580 and p+ 2.580. 


Again, these results hold whatever the values of the mean : and the 
standard deviation o. The results are illustrated in Figure 6.6. 


p-1.640 pe p+ 1.64o 1960 p p+1.960 p-2580 pp + 2.580 
Figure 6.6 Areas under normal curves 


These general results are of practical importance and, as you will see, we 
shall make use of the second result in particular in the next two chapters. 


Example 6.2 Heights of Cambridge men 
According to the normal distribution used to model the variation in the 


heights of Cambridge men in 1902, within what ranges were the heights of 
approximately 90%, 95% and 99% of Cambridge men? 
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Solution 


Applying the general results above to the normal model for the heights of 
Cambridge men in 1902, we have the following: 
© approximately 90% of Cambridge men were between 
68.9 — 1.64 x 2.57 ~ 64.7 and 68.9 + 1.64 x 2.57 ~ 73.1 inches tall; 
© approximately 95% of Cambridge men were between 
68.9 — 1.96 x 2.57 ~ 63.9 and 68.9 + 1.96 x 2.57 ~ 73.9 inches tall; 


© approximately 99% of Cambridge men were between 
68.9 — 2.58 x 2.57 ~ 62.3 and 68.9 + 2.58 x 2.57 ~ 75.5 inches tall. 


According to the normal model used in Activity 6.1, within what ranges 
are the lengths of the following percentages of all cuckoo eggs? 


(a) 90%  (b) 95% = (c) 99% 
Comment 


A solution is given on page 46. 


All the results about the proportion of values within a given number of 
standard deviations of the mean are special cases of the following general 
result, 


If a normal distribution is used to model the variation in a 
population, then, according to the model, the proportion of the 
population within k standard deviations of the mean is the same 
whatever the values of the mean jz and the standard deviation o. 


This result means that the area shaded in Figure 6.7 depends on the value 
of k, but not on the values of js and a. (See the box on page 39 for a 
suggestion on how you could check this result.) 


u-ke p pt+ko 
Figure 6.7 k standard deviations from the mean 


In particular, for any values of and a, this area is equal to the area 
between —k and k under the curve of a normal distribution with mean 
zero and standard deviation one (j= 0 and o = 1). This is illustrated in 
Figure 6.8. 
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Figure 6.8 Equal areas 


SECTION 6 PROPERTIES OF NORMAL DISTRIBUTIONS 


The normal distribution with mean equal to zero and standard deviation 
equal to one is called the standard normal distribution. A sketch of 
the standard normal curve is shown in Figure 6.9. 


af HBr a2” 88 
Figure 6.9 The standard normal distribution 


Tables giving numerical areas under the standard normal curve are 
included in many statistics textbooks and in books of statistical tables. 
Since the area within k standard deviations of the mean is the same for all 
normal distributions, these tables can be used to find areas under any 
normal curve. Even though statistics software packages are available to do 
calculations for normal distributions, printed tables are still used when it is 
not convenient, or not worthwhile, to switch on a computer. If you study 
statistics in the future, you will almost certainly learn how to use such 
tables. However, in MST121, the computer will be used for all calculations 
involving normal distributions. 


Checking a result 

If you are interested, you could check the result concerning the 
proportion of values that are within & standard deviations of the 
mean in a normal distribution by finding the area under a normal 
curve between ps — ko and + ko. This area is given by the integral 


1 wether Ne 
exp = (- a) dx. 
oV lt Junko 2 o 
You can use Mathcad to evaluate this integral. You should find that 
the result given by Mathcad contains k but does not involve j and o, 


Summary of Section 6 


In this section, some properties of normal distributions have been 
summarised. The main result that will be used in the next two chapters is 
as follows. 

For a population modelled by a normal distribution, the proportion 

of the population within k standard deviations of the mean (that is, 

between ps — ko and ps + ko) is the same whatever the values of 

and o. 
In particular, approximately 95% of the population are within 1.96 
standard deviations of the mean, that is, between ps — 1.960 and ps + 1.960. 


You do not need OUStats for 
this exercise. 
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CHAPTER D2 MODELLING VARIATION 


Exercise for Section 6 
Exercise 6.1 Sons’ heights 


In Section 4, you explored Pearson’s data on the heights of father-son 
pairs. For fathers who were 67 inches tall, the mean height of their sons 
was 68.0 inches and the standard deviation was 2.21 inches. Suppose that 
a normal distribution with parameters 4 = 68.0 and o = 2.21 is used to 
model the heights in inches of the sons of fathers who are 67 inches tall. 


(a) Calculate a range of values within which, according to the model, the 
heights of 95% of sons of 67-inch-tall fathers lie. 


(b) A 67-inch-tall man has a son who is 63 inches tall. How many 
standard deviations below the mean height for sons of 67-inch-tall 
fathers is this son, according to the model? Another 67-inch-tall man 
has a son who is 72 inches tall. How many standard deviations above 
the mean height for sons of 67-inch-tall fathers is this son? Would you 
regard either of these sons as unusually short or tall for the sons of 
67-inch-tall fathers? 


Summary of Chapter D2 


In this chapter, the idea of using a curve to model the variation observed 
in data has been introduced. One particular model — the normal 
distribution — has been discussed in some detail. But you should remember 
that this is not the only possible model. In Activity 1.2, you saw examples 
where differently-shaped curves are needed. There are other distributions 
which are widely used for modelling variation, and if you study statistics in 
the future you will meet some of them. However, the normal distribution is 
of great practical importance, as you will see in the next chapter. 


You have also seen that sample statistics, such as the sample mean 7 and 
the sample standard deviation s, may be used to estimate population 
parameters, such as the population mean j and the population standard 
deviation a. A large part of your work has involved learning to use the 
data analysis software OUStats, and you have been encouraged to apply 
the skills you have acquired using Mathcad to aid you in this. 


In the next two chapters, you will be applying the idea of using sample 
statistics to estimate population parameters; you will also be using some 
of the properties of the normal distribution discussed in this chapter. 
Additionally, some further features of OUStats will be introduced. So, 
before you move on to the next chapter, it is important for you to reflect 
on what you have learned. To help you focus your attention on some of the 
important features of this chapter, try the activities suggested below. 


Activity 7.1 Modelling variation: a presentation 


Suppose that you have been asked to give a short presentation to show 
what is meant by the statement: ‘the variation in the weights of week-old 
nestlings may be modelled by a normal distribution’. 


Prepare a short presentation (10-15 minutes) to show how such a model 
could be used to calculate probabilities or answer questions about the 
proportions of nestlings of different weights. 


(You will find some hints on preparing presentations on the reverse side of 
Learning File Sheet 3.) 


Activity 7.2 Learning to use OUStats 


In this chapter, you have been introduced to the main features of OUStats, 
the data analysis part of the statistics software. If you have mastered the 
facilities used in this chapter, then you should find using the remaining 
ones straightforward. How easy or difficult did you find it to learn to use a 
new Windows-based software package? Did your experience using Mathcad 
help you with OUStats and, if so, in what way did it help you? You might 
find it helpful to check the notes you made at the end of Section 3 on the 
similarities and differences between the two packages. Is there anything 
you wish to add to your notes as a result of your work on Sections 4 and 5? 
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CHAPTER D2 MODELLING VARIATION: 


Learning outcomes 


You have been working towards the following learning outcomes. 


Terms to know and use 
Mode, unimodal and bimodal, skewed, left-skewed and right-skewed, 
probability density function, normal curve, normal distribution, 
sample mean and sample standard deviation, population mean and 
population standard deviation, sample statistics and population 
parameters, 


Symbols and notation to know and use 
The notation Z for the sample mean and s for the sample standard 
deviation. 
The symbol x for the mean of a probability distribution or for a 
population mean, 
The symbol ¢ for the standard deviation of a probability distribution 
or for a population standard deviation. 


Ideas to be aware of 


© How the uncertainty about the value of a continuous variable can be 
represented by a curve, that is, by a probability density function. 

© How to interpret an area under the graph of a probability density 
function (such as that of the normal distribution) as a probability or 
as a proportion. 

© How the sample statistics Z and s may be used to estimate the values 
of the population parameters 1 and o. 


Features of OUStats 


© Open and close an existing data file, and obtain information about the 
data; obtain summary statistics for the data in a file. 


© Obtain a frequency diagram for a sample of data, and fit a normal 
curve to the data. 


© Find areas under a normal curve. 

© Generate random samples from a normal distribution. 
© Print output. 

© Obtain a scatterplot for paired data. 

Learning skills 


© Observe similarities between software packages, and apply skills 
acquired with one package to help to learn another. 


Investigating processes to aid understanding 
© Conjecture a general result from results obtained in special cases. 


Appendix: Integrals with infinite limits and the 
normal distribution 


In Section 1, the area between a normal curve and the z-axis is defined as 
an integral with infinite limits. And in Section 2, the mean and variance of 
a continuous distribution with probability density function f are defined as 
integrals with infinite limits. This appendix explains briefly one approach 
to defining an integral with infinite limits, and looks again at the integrals 
introduced in Sections 1 and 2. 


1. Integrals with infinite limits 
Suppose that, for any number N > 0, the definite integral 


[seer 


can be calculated. If the value of this integral approaches a limiting value 
as N becomes larger and larger, then the expression 


[sear 


is defined to be equal to this limiting value; that is, 
oo N 
i. f(x) de = lim j f(a) de. 
= N00 Jy 


2. The total area under a normal curve 


For a normal distribution with probability density function f given by 


f(z) = ae (-} (5*)) (-00 < @ < 00), 


the total area between the normal curve and the w-axis is given by 


[sae 


If you are interested (and have the time), you could investigate this integral 
using Mathcad, as follows. For any y and @ of your choice, evaluate 


ale (-} (5*)) dx 


for successively larger values of N. 
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CHAPTER D2 MODELLING VARIATION 


3. The mean and standard deviation of a normal distribution 

In Subsection 2.1, it was stated that the mean of a normal distribution is 
given by the value of the parameter jz. This result can be deduced using 

the symmetry of the normal curve, but you might like to use Mathcad to 
investigate the integral 


[ 282) de= et ad (-} (4)) dr. 


In Subsection 2.2, it was stated that the standard deviation of a normal 
distribution is given by the parameter a. The variance is defined by the 
integral 


[ije-mieren fete (3 (254)) an 


and the standard deviation is the square root of the variance. You could 
also use Mathcad to investigate the above integral, keeping ps and o as 
constants. 


Solutions to Activities 


Solution 1.2 


For the earthquake data in Figure 1.4(a), short times 
are the most common, and longer times occur less 
often. The frequency diagram has a long right tail. 

A sketch of a possible model for times between major 
earthquakes is given in Figure S.1. 
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Figure S.1 A possible model for times between 
earthquakes 


The frequency diagram for the wing lengths of 
first-year meadow pipits has a peak at about 81mm, 
and there is a suggestion that there might be a 
second smaller peak at about 76 or 77mm. 

A possible model for wing length is sketched in 
Figure S.2. 
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Figure S.2 A possible model for wing lengths of 
meadow pipits 


The frequency diagram for the durations of eruptions 
of Old Faithful geyser has two peaks or modes, at 
approximately 1.8 and 4.0 minutes; it is bimodal, 

A sketch of a possible model is shown in Figure 8.3. 
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Figure S.3 A possible model for the durations of 
eruptions of Old Faithful geyser 


The frequency diagram for the women’s weekly 
earnings data has a single peak at about £200. It has 
a longer right tail than left tail. A possible model for 
the variation in weekly earnings is shown in 

Figure $.4. 
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Figure S.4_ A possible model for women’s weekly 
earnings in April 1994 


Solution 1.4 


(a) The shaded area represents the proportion of all 
Cambridge men in 1902 who were between 69 
and 71 inches tall. Alternatively, it represents 
the probability that a man selected at random 
from all Cambridge men would have been 
between 69 and 71 inches tall. 


(b) The shaded area represents the proportion of 
all Cambridge men in 1902 who were under 
65 inches tall, Alternatively, it represents the 
probability that a man selected at random from 
all Cambridge men would have been less than 
65 inches tall. 


(c) The shaded area represents the proportion of all 
Cambridge men in 1902 who were over six feet 
tall. Alternatively, it represents the probability 
that a man selected at random from all 
Cambridge men would have been over six feet 
tall. 


Solution 6.1 


According to the model, the lengths of about 99.7% 
of all cuckoo eggs should be between js — 30 and 

4. + 30; that is, between 22.4 —3 x 1.08 ~ 19.2mm 
and 22.4 +3 x 1.08 ~ 25.6mm. 


For the sample of 243 eggs, no egg was as long as 
25.6 mm and no egg was shorter than 18.75 mm. 
One was recorded as 19mm long (between 18.75 mm 
and 19.25mm), so the length of at most one of the 
243 eggs was outside the range 19.2mm to 25.6 mm. 
So ‘almost all’ the eggs were within three standard 
deviations of the mean. 


Solution 6.2 


(a) The lengths of 90% of eggs should lie between 
22.4 — 1,64 x 1.08 ~ 20.6mm and 
22.4 + 1,64 x 1.08 ~ 24.2mm. 


(b) The lengths of 95% of eggs should lie between 
22.4 — 1.96 x 1.08 ~ 20.3mm and 
22.4 + 1.96 x 1.08 ~ 24.5mm. 


(c) The lengths of 99% of eggs should lie between 
22.4 — 2.58 x 1.08 ~ 19.6 mm and 
22.4 + 2.58 x 1.08 ~ 25.2mm. 


Solutions to Exercises 


Solution 2.1 
(a) The sample statistics are 


E=3, s~ 2.345 20788 ~ 2.3, 


The details of the calculations are given below 
for you to check your working if your calculator 
does not have statistical facilities for the mean 
and standard deviation. 


2 2-F (x-)? 
2 =1 ie 
ig 4 16 
3 0 0 
1 =2 4 
Cost 1 
15 22 
From the first column of the table, we can 
calculate 
E=15/5=3. 


Then the values in the second and third columns 
can be calculated, to give 


8 = 1/22/4 ~ 2.3, 


(b) The sample mean % may be used to estimate the 
mean number of items of mail delivered to each 
house on the estate: 3. So an estimate of the 
total number of items of mail delivered to the 
estate on that morning is 3 x 3000 = 9000, 


Solution 2.2 
The sample statistics are 
E=507, s =~ 157.297171 ~ 157.3, 


The details of the calculations are given below. 


x z-% (#-z)? 


310 —197 38809 
635 128 = 16384 
464 —43 1849 
520 13 169 
381-126 15876 
732 225 50625 
3042 123712 


From the first column of the table, we can calculate 
% = 3042/6 = 507. 


Then the values in the second and third columns can 
be calculated, to give 


8 = /123712/5 ~ 157.297 171 ~ 157.3. 


So the mean gross weekly earnings in 1995 of the 
sample of six mechanical engineers was £507, and 
the sample standard deviation of gross weekly 
earnings was approximately £157.30. 


Solution 6.1 


(a) According to the model, the heights of 95% of 
sons will be within 1.96 standard deviations of 
the mean, that is, between 
68.0 — 1.96 x 2.21 ~ 63.7 inches and 
68.0 + 1.96 x 2.21 ~ 72.3 inches. 


(b) At 63 inches tall, the first son is 5 inches shorter 
than the mean height (68 inches) of sons of 
67-inch-tall fathers, that is, 5/2.21 ~ 2.26 
standard deviations below the mean. 


At 72 inches tall, the second son is 4 inches 
taller than the mean height of sons of 
67-inch-tall fathers, that is, 4/2,21 ~ 1.81 
standard deviations above the mean. 


One way to look at the last question is as 
follows. Since the heights of 95% of sons are 
within 1.96 standard deviations of the mean 
height, and only 5% are 1.96 or more standard 
deviations from the mean height, we might 
regard a son whose height is at least 1.96 
standard deviations above the mean as unusually 
tall and a son whose height is at least 1.96 
standard deviations below the mean as unusually 
short. On this basis, the first son is unusually 
short; but we would not consider the second son 
to be unusually tall, as his height is less than 
1.96 standard deviations above the mean. 
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